Glass transition as a decoupling-coupling mechanism of rotations 
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We introduce a three-dimensional lattice gas model to study the glass transition. In this model 
the interactions come from the excluded volume and particles have five arms with an asymmetrical 
shape, which results in geometric frustration that inhibits full packing. Each particle has two degrees 
of freedom, the position and the orientation of the particle. We find a second order phase transition 
at a density p « 0.305, this transition decouples the orientation of the particles which can rotate 
without interaction in this degree of freedom until p — 0.5 is reached. Both the inverse diffusivity and 
the relaxation time follow a power law behavior for densities p < 0.5. The crystallization at p = 0.5 
is avoided because frustration lets to the system to reach higher densities, then the divergencies are 
overcome. For p > 0.5 the orientations of the particles are coupled and the dynamics is governed 
by both degrees of freedom. 

PACS numbers: 64.70.Pf 
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I. INTRODUCTION 

In the last years, a great deal of work has been done to 
obtain a fundamental understanding of the glass transi- 
tion. Many questions about the equilibrium and the dy- 
namical properties of the glassy state remain ananswered. 
It is not clear if there is a true phase transition and what 
is the role that geometric frustration plays on it. The re- 
lations between the equilibrium and the dynamical prop- 
erties are not understood Q| . The mode coupling theory 
for supercooled liquids || predicts the existence of a tem- 
perature T c at which there is a crossover from a liquid to 
a glassy state. In the glassy state the dynamics would be 
dominated by complex activated processes. For tempera- 
tures T > T c but close to T c there is a power law behavior 
of the relaxation time and also of the inverse diffusivity. 
It is not understood whether T c is a purely kinetic transi- 
tion temperature || or if it is a true thermodynamic glass 
transition which is kinetically avoided || . Several lattice 
gas models have been used to simulate glassy systems and 
have reproduced some aspects of the glassy phenomenol- 
ogy. For example the Hard Square Model (HSM) j| , the 
Kinetically Constrained Model || [| , the Frustrated Ising 
Lattice Gas Model [0 , the one introduced by Ciamarra et 
al. ||, and recently the Lattice Glass Model (LGM) ||. 
The LGM model with density constraint / = is equiva- 
lent to the the three-dimensional HSM model also called 
Hard Cubic Model (HCM). The LGM model relates the 
glass transition to a first order phase transition. 

In this paper we consider a three-dimensional lattice 
gas model, which contains as main ingredients only geo- 
metric frustration without quenched disorder and with- 
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FIG. 1: Schematic picture of one particular configuration in 
a system size of N = 2 3 and density p = 0.5. Shadow spheres 
represent holes in the system. Black spheres are particles with 
five arms each one. 



out kinetic or density constraints, as quenched disorder 
is not appropriate to study structural glasses and kinetic 
or density constraints are some how artificial. Similar 
models have already been proposed and studied in two- 
dimensional systems [ll], |l2| and applied to study 
granular material p3|. 



II. THE MODEL 

Our model is a generalization to three dimensions of 
the two-dimensional model studied in Ref . [Jl2j . It can be 
considered as an illustration of the concept of frustration 
arising as a packing problem. We have particles with five 
arms and they occupy the vertices of a cubic lattice with 
one of six possible orientations. Assuming that the arms 
cannot overlap due to excluded volume, we see that only 



for some relative orientations two particles can occupy 
nearest-neighbor vertices. Consequently, depending on 
the local arrangement of particles, there are sites on the 
lattice that cannot be occupied (see Fig. [l]). This type of 
"packing" frustration thus induces defects or holes in the 
system. We impose periodic boundary conditions in the 
cubic lattice of size TV = L 3 . The maximum of density is 
Pmax = 3/5 at which all possible bonds are occupied by 
an arm. Here we have two degrees of freedom for each 
particle, the position and the orientation of the particle. 
This model is the HCM model when the particles have 
six instead five arms. Our model would be also similar 
to the LGM model with the density constraint 1 = 1, Wc 
will compare the results found in our model with the ones 
obtained in these two models. We have used two algo- 
rithms in order to make the simulations. The first one 
(CA) is the Monte Carlo simulation at fixed density in 
the canonical ensemble, we have simulated the diffusion 
and rotation dynamics of the particles by the following 
algorithm: i) Pick up a particle at random; ii) Pick up 
a site at random between the six nearest neighbor ones; 
iii) Choose randomly an orientation of the particle; iv) If 
it does not cause the overlapping of two arms, move the 
particle in the given site with the given orientation; v) 
If the diffusion movement is not possible, choose a ran- 
dom orientation and try to rotate the particle to this new 
orientation; vii) Advance the clock by 1/N, where TV is 
the number of sites, and go to i). The second algorithm 
(GCA) is the grand canonical ensemble, the diffusion and 
rotation dynamics is as in the CA simulations but now 
a reservoir with chemical potential p is coupled to each 
lattice site which can create (if it does not cause the over- 
lapping of two arms) or destroy particles. As we expect 
the GCA simulation reaches the equilibrium faster. We 
will use the GCA simulation in order to find the behavior 
of the density p with p || . 



III. RESULTS 
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FIG. 2: Diffusion constant D as a function of density p for a 
system size of TV = 14 3 . The fitting function is a power law 
D = 2.11(0.523 - p) 2 45 . Inset: The inverse of the density 
as function of p on a lattice of size TV = 14 3 . The data arc 
obtained from the GCA with 2xl0 7 Monte Carlo Steps per 
Particle at a fixed increasing rate of the chemical potential 
between p = and p = 27. 

7 = 2.45 ± 0.01 (see Fig. ||). This anomalous behavior 
of D near of p c would indicate a crossover density, from 
liquid to glass phase, where activated processes dominate 
in the glass phase. As it happens in the HSM model 
the finite size effects for the diffusivity in the HCM model 
are very large when we do CA simulations, it is because 
the particles can be enclosed in cages and the diffusion is 
blocked. It prevents to reach densities close to the max- 
imum density. Nevertheless, in our model the finite size 
effects are only important for p > 0.58 when TV = 14 3 . 
This is because our model has two degrees of freedom 
and rotations prevent to find blocked configurations for 
p < 0.58. 



We first study a possible first order phase transition 
in our model. In the inset of Fig. the inverse of the 
density is plotted as a function of the chemical poten- 
tial. As in Ref. |J we make GCA simulations to obtain 
this figure. A maximum of density very close to p max is 
reached without any discontinuity, although we observe 
finite-size effects which prevent to reach p max for the lat- 
tice sizes studied here. So, first order phase transition is 
not present in our system. In the HCM model we also 
observe similar behavior of the density with p reaching 
the maximum of density continously, p max = 0.5 in the 
HCM model. Instead, in the LGM model with I > 1 
there is always a first order phase transition 0. 

We now calculate the diffusion coefficient D from the 
mean-square displacement of the particles at very long 
times with the CA simulations. The values obtained for 
D are well fitted by a power law close to p c = 0.52±0.005 
and for densities lower than p = 0.5, D oc (p c — p) 1 with 



In order to understand what happens at p f=a 0.5 we 
study the following microscopic order parameter. As in 
anti-ferromagnetic systems, the cubic lattice is divided 
into two interpenetrating sublattices (A and B) , a site in 
a sublattice has six nearest neighbor sites which belong 
to the other sublattice. The order parameter is defined 



PA - Pb 
Pa + Pb 



(1) 



where pa and pb are the equilibrium concentrations of 
the particles in the sublattices A and B and we have 
P = PA + Pb- This parameter can be used to study con- 
centration sublattice ordering. When the particles prefer 
to stay in one of these sublattices then <j> ^ 0. In the left 
inset of Fig. |^ we show <j> as a function of the density. 
We can see that it is different to zero for p > pf w 0.3 
and it increases until p « 0.5, then there is a maximum. 



FIG. 3: Concentration of particles in sublattices A, pA (as- 
terisks), and B, pb (solid circles), as a function of the density 
p for a system size of L — 14. Left inset: Order parameter <f> 
as a function of the density p. Right inset: Compressibility « 
as a function of the density p. 



For higher densities it decreases linearly with the density. 
In Fig. ^ we show the concentration in both sublattices. 
We see that at pi the concentrations pa and ps begin to 
be different each other, the particles prefer to stay in a 
sublattice. The concentration in the sublattice A has the 
maximum value when p = 0.5 is reached. A sublattice 
is full of particles at this concentration while the other 
one is in practice empty. The particles begin to occupy 
the empty sublattice for densities higher than p = 0.5 
remaining the other sublattice full, then the parameter 
4> decreases. Here we observe that there are not frozen 
particles at these densities, pa and ps are equilibrium 
concentrations. The order parameter <p can take positive 
and negative values, it depends on which sublattice has 
higher density for p > pf, in Fig. [| we have tj> > for 
p > pi because pa > Pb- Similar behavior for the pa- 
rameter (f> is found in the HCM model but pi « 0.22 is 
lower than the one obtained in our model and the system 
crystallizes at p — 0.5, its maximum of density, then the 
particles are frozen in a sublattice and it is not possible 
to reach higher densities. So, the last part of Fig. || (for 
p > 0.5) is not found in the HCM model. As we will sec 
below at pf there is a continuous phase transition in our 
model and also in the HCM model. In the LGM model 
with 1 = 1, we have found a first order phase transition 
with a discontinuity of the density, when is plotted as 
a function of p, from a density p w 0.32 to a density 
p = 0.5, then the system crystallizes and the particles 
are frozen in a sublattice. 



We now study the equilibrium second order phase tran- 
sition in our system. For that, we define the compress- 



FIG. 4: Finite size scaling of \ an d 5 (inset) for lattice sizes 
L = 8 (•), 10 (■), 12 (♦), 14 (A), and 16 {<). p f = 0.305 ± 
0.005, r) = 0.04 ± 0.01, and v = 0.63 ± 0.05. 



ibility by the following expression 

*=4£«^>-<^> 2 ) 



(2) 



where rij = 0, 1 is the occupation number of site i and 
(• • •) indicates equilibrium average. In the right inset of 
Fig. U we see the behavior of k with the density. We 
find that the compressibility is different to zero for den- 
sities higher than pi w 0.305 and it has a maximum at 
p = 0.5, decreasing for higher densities. This behavior 
is similar to the one found for the parameter ip. The as- 
sociated susceptibility \ is given by \ = N ((ft 2 ) — (k) 2 ) 
and the Binder's cumulant is g = (3 — (k 4 )/(k 2 ) 2 )/2. 
Around a continuous phase transition \ an( l .9 should 
obey the finite-size scaling x{p) = L 2 ~ v x \L 1 ^ u {p — Pf)l 
and g(p) = g \L x l v {p — pf)] where x [ x ] an d 9 N are uni- 
versal functions and pi is the critical density. From the 
finite size scaling (see Fig. ^) we find a continuous phase 
transition at pf = 0.305 ± 0.005 which belongs to the 
three-dimensional Ising universality class, r\ = 0.04±0.01 
and v = 0.63 ± 0.05. In the HCM model we have found 
a second order phase transition which also belong to the 
same universality class but with pf « 0.22. In Fig. || we 
can see that the diffusion constant is not affected by the 
second order phase transition. 

We now study the relaxation of the autocorrelation 
functions of the density fluctuations 



(lp q l 2 > 



(3) 



where (■ • •) denotes average over the reference time t' and 





FIG. 5: Correlation functions of the density fluctuations 4 > q (t) 
for q = (n/3, 7r/3, 7r/3), system size L — 14, and densities 
(from bottom to top) p = 0.45, 0.46, 0.47, 0.48, 0.49, 0.5, 
0.505, 0.51, 0.515, 0.52, 0.53, and 0.54. The solid lines are 
fitting functions corresponding to stretched exponential func- 
tions. For p > 0.515 we have /3 = 0.94 ± 0.01. Inset: D as a 
function of the relaxation time r. The arrow shows the density 
p — 0.5. The solid line is the fitting function D = 0.135/t. 

p q is the Fourier transform on the lattice of the density 

1 ™ 

where ti(t) is the position of the ith particle at time t, 
n is the number of particles and q is the wave number. 
Because of the periodic conditions on the cubic lattice 
q = (27T / 'L)(n x , n y , n z ), with n x , n y , n z = 1, • • • , L/2. 
Fig. | shows $ q (i) for q = (tt/3, tt/3, tt/3) and differ- 
ent densities. We can see a two-step relaxation decay for 
p > 0.5, the second relaxation step can be fitted by a 
stretched exponential form, f(t) = /oexp [(t/r) 13 ] where 
the exponent (3 = 0.94 ±0.01 remains constant with p for 
p > 0.515. For p < 0.49 we only have a one-step relax- 
ation decay. The relaxation time r can be obtained from 
$ q (t). We find that it is proportional to the inverse of 
the diffusivity, r oc D , for the whole range of densities 
studied here (see inset of Fig ||). Thus, the relaxation 
time follows a power law at p < p c with the same expo- 
nent than the one obtained in the power law of D _1 (see 
Fig. §. 

We now present the results for the self-part of the au- 
tocorrelation function of the density fluctuations, defined 

as 

n« = ^E< e * q(r!(t ' +t) ~ ri(0) >' (5) 

i 

where (• • •) denotes average over the reference time t' . 
Fig. H shows $ q (t) corresponding to q = (tt, 0, 0) for 
densities p = 0.35, 0.45, 0.5, and 0.53. For the whole 



FIG. 6: Self-part 3>q(t) of the density-density autocorrelations 
for q = (-7T, 0, 0) and densities p = 0.35, 0.45, 0.5, and 0.53. 
Dotted lines are fitting functions corresponding to stretched 
exponential functions. Inset: Parameter /3 of the stretched 
exponential function for $q(£) as a function of the density 
(for each density we plot the error bar of (5). 

range of densities studied here we find that the whole 
time interval of $ q (i) can be fitted by a stretched expo- 
nential function where the exponent (3 depends on the 
density. In the inset of Fig. |] we show (3 as a function of 
the density. The exponent (3 decreases with the density 
until a density near pf is reached. Starting from this den- 
sity, which corresponds to the second order phase tran- 
sition, the exponent increases until p — 0.5 is reached. 
For p > p c it becomes constant (within the error bars) 
(3 « 0.93. The relaxation time obtained from the fit of 
$ q (t) is proportional to the inverse of the diffusivity. 

In order to study the role of the orientation of the par- 
ticles we define a self-overlap parameter similar to that 
defined in jl^] for liquids but which also takes into ac- 
count the orientation of the particles, besides their posi- 
tion. The orientation of a particle is defined by the dis- 
crete values of the two orthogonal angles Qi = 0, tt/2, tt, 
or 37r/2 and <fi = 0, tt/2, tt, or 37r/2. We define the 
self-overlap as 

1 n 

<?(*) = -jvrX!fa<(*')"i(* ' + t) 



cos [9i(t' +t)- 90)} cos [ Vi {t' +t) - <pi (*')]) , (6) 

here (■ • ■) denotes average over the reference time t'. If all 
the particles have theposition and the orientation frozen 
then q(t) = 1. Fig. shows the parameter q(t) for dif- 
ferent values of the density. The plateau becomes visible 
for densities higher than p w 0.5. From this density there 
are an important number of particles which have frozen 
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FIG. 7: Relaxation functions of the self-overlap q(t) for sys- 
tem size L — 14 and densities (from bottom to top) p = 0.45, 
0.46, 0.47, 0.48, 0.49, 0.5, 0.51, 0.515, 0.52, 0.53, 0.54, 0.55, 
0.56, 0.57, and 0.58. 

the position and the orientation for a long time. The 
number of frozen particles and the time during they are 
frozen increase with the density. We can fit the second 
relaxation step with the stretched exponential function, 
but now the exponent (3 decreases with the density from 
f3 = 0.94 for p = 0.52 until (3 = 0.77 for p = 0.58. ' 



rotate freely. For densities higher than p = 0.5 geomet- 
ric frustration begins to work and rotations are governed 
by complex collective processes. Then, the two degrees 
of freedom arc important in the diffusivity movement of 
the particles. For p < 0.5 the system is going to a crys- 
talline state with all the particles frozen in a sublattice, 
this brings to a power law divergency of the relaxation 
time and the inverse of diffusivity for p < 0.5. But frus- 
tration lets to the system reach higher densities and crys- 
tallization is avoided and the divergencies are overcome. 
Then, vibrational effects are observed which bring to the 
two-step relaxation decay in the density correlations and 
in the self-overlap parameter. Thus, the glass transition 
is purely a kinetic transition in our model. Geometric 
frustration plays a fundamental role, without frustration 
the arrest would be close to p c , but also the second or- 
der phase transition is very important, which decouples 
the orientation of the particles. In the two-dimensional 
model which has not second order phase transition, 
we do not observe these anomalies in the diffusivity and 
relaxation time. The order parameters <fi and k exhibit a 
maximum at the glass transition p w 0.5. 

We have found that the diffusion constant is not af- 
fected by the second order phase transition. The self-part 
of the autocorrelation function of the density fluctuation 
can be fitted by a stretched exponential function with 
an exponent (3 that has a minimum value at the second 
order phase transition and a local maximum at the glass 
transition. 



IV. CONCLUSION 

We have proposed a three-dimensional lattice gas 
model, based on the concept of geometric frustration 
which is generated by the particle shape. In this model 
a second order phase transition decouples the orienta- 
tion of the particles which can rotate without interac- 
tion in the orientation degree of freedom until p = 0.5 is 
reached. This is because in practice the particles remain 
all the time in a sublattice and then the particles can 
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